Continuous Observation of Interference Fringes from Bose Condensates 



J. I. Cirac, 1 C.W. Gardiner, 2 M. Naraschewski, 3,4 and P. Zoller 5 
1 Departamento de Fisica Aplicada, Universidad de Gastilla-La Mancha, 13071 Giudad Real, Spain 

2 Physics Department, Victoria University, Wellington, New Zealand 
3 Max-Planck-Institut fur Quantenoptik, Hans-Kopfermann-Strafie 1, D-85748 Garching, Germany 
4 Sektion Physik, Universitat Munchen, Theresienstrafie 37, D-80333 Munchen, Germany 
5 Institut fur Theoretische Physik, Universitat Innsbruck, Technikerstrafie 25, A- 6020 Innsbruck, Austria 

(June 3, 1996) 



We use continuous measurement theory to describe the 
evolution of two Bose condensates in an interference experi- 
ment. It is shown how in a single run the system evolves into 
a state with a fixed relative phase, without violating particle 
number conservation. 
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Recent observation of Bose-Einstein Condensation 
(BEC) [1-3] has initiated theoretical discussions regard- 
ing the properties of Bose condensates [4]. Of particular 
interest have been questions related to the phase of the 
condensate [5] . The assumption of such a phase as a re- 
sult of a broken gauge symmetry allows for a natural ex- 
planation of many physical phenomena, but also implies 
that the state of the condensate is a linear superposi- 
tion of states with different particle number. However, 
in second quantized formalism of nonrelativistic Quan- 
tum Mechanics, as is usually employed to describe BEC, 
all observables commute with the atomic number opera- 
tor N which thus plays the role of a superselection rule. 
As a consequence, starting from a state with fixed atomic 
number or a mixed state, which is diagonal in the atomic 
number basis, no atomic superpositions and coherences 
will develop [6] . 

In order to resolve this seeming contradiction, we use in 
this paper the language of continuous measurement the- 
ory [7] to describe a single realization of an interference 
experiment between two independent condensates. We 
will discuss how the state of the two condensates evolves 
as atoms are detected. In particular, from our analy- 
sis it follows how a state of well defined relative phase 
builds up dynamically in a single experimental run as a 
consequence of the von Neuman projection postulate of 
quantum theory. We emphasize that our description does 
not contradict atomic number superselection rules, and 
allows one to understand the coexistence of both particle 
number conservation and the phase of a condensate in a 
general situation. The problem addressed in the present 
paper is related to very recent work by Javanainen and 
Yoo [8], and Naraschewski et al. [9], where it is shown 
that two independent Bose condensates prepared in Fock 
states may form a measurable interference pattern. How- 
ever, here we will show how the interference pattern is 
formed dynamically, and, what is more important, how 
the state of the two condensates collapses during the pro- 



cess of detection of the many individual atoms 

condensate 1 f~ 

c V7 



condensate 2 



• _ • si 



scree 




FIG. 1. Schematic setup of interference of two Bose conden- 
sates. 



We consider the situation depicted in Fig. 1: two sta- 
tistically independent Bose condensates A and B of iden- 
tical atoms are dropped and detected by atomic counters 
in the z = plane. For simplicity, we assume that all 
the particles in the first (second) condensate occupy the 
same mode, which is a plane wave of wavevector k a (k b ). 



We also assume that k 
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they have opposite momenta in the x direction [10]. The 
initial density operator describing the state of the con- 
densates can be written as po = p^ <X> p^ . Consistent 
with the atomic number superselection rule we require 
to be diagonal in the atomic particle numbers, 



p ij) = ^2ptt ) \n,j)(n,j\ (j = a,b) 



This includes as a special case Fock states p^' = 
\n, j)(n, j\, i.e. states with a fixed atom number (for 
the existence of interference fringes in single realization 
of experiments with Fock states see Refs. [8,9]). The 
derivations which follow are considerably simplified [11] 
if we assume for the initial states of the condensates an 
ensemble of phase averaged coherent states with fixed 
amplitudes a and b, respectively, 
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which correspond to Poissonian atomic number distribu- 
tions pn ^ = e~i j 2n /n\, (j = a,b). 

As the atoms are detected, the number of atoms in each 
condensate decreases. Due to the fact that the atoms are 
identical, when one atom is detected one cannot deter- 
mine whether it belonged to the first or the second con- 
densate, which gives rise to interference. In other words, 
the number of atoms in each of the condensates is not 
conserved. However, the total number of atoms must be 
conserved, and thus the density operator describing the 
whole system must remain diagonal in the total number. 
It is therefore convenient to define 

R ab (<f>) = P ^t±\ae^){ae^\® \be^-^){be^~^\ 
Jo 27T 

(3) 

as the density operator of a state with fixed relative phase 
<f>. Note that this state is diagonal in the total atom 
number (consistent with overall number conservation), 
but is an entangled state with respect to the particles in 
both condensates [12]. The initial state (2) can be then 
rewritten as 

po = r tt R ^ = r ^w^w w 

Jo 2lr Jo 2lr 

which corresponds to a state with uniformly distributed 

random relative phase (f^J (4>) = 1). 

In order to describe the evolution of the density op- 
erator as the atoms are detected at the z = plane we 
use a simple continuous measurement model based on the 
master equation [13] 

p = —iuj[a) a + b'b, p] + n(2apa^ — a) dp — pa^ a) (5) 
+n{2bptf - Pip- pb r b), 

where a and b are annihilation operators for particles in 
the first and second condensate, respectively. The evolu- 
tion frequency u> and the loss rate due to detection k are 
assumed to be equal for the two modes, in accordance 
with their corresponding momenta k a and k b . Note that, 
according to this master equation, the density operator 
does not develop atomic superpositions nor coherences. 
The relationship between this master equation and the 
physical situation we are considering becomes apparent 
if we rewrite (5) as 

r-2ir 

p = -iH eff p + ipHlff + / d(f)J^p, (6) 
Jo 

where 

J r pp = 2 K (d + be^)p(a + be i ^, ^e[0,27r), (7a) 
H eff = (u - i K )(a t a + & t &), (7b) 

are the (nonhermitian) effective Hamiltonian, and recy- 
cling superoperator, respectively. According to contin- 
uous measurement theory [7], the interpretation of (6) 



in the present situations is as follows: a detection of an 
atom at the position xo is associated with the action of 
the jump superoperator as given in (7a), 

Pt N +st oc J^p tN (8) 

with <f> = [k b x — k£]xo- That is, an atom is removed from 
the first or the second condensate, and the position of 
the atom is associated with a relative phase <f> of the 
plane waves. The time evolution between two subsequent 
detections is governed by 

Pt N+1 oc e - iH MtM +1 -t N ) ptNe iHl ff (t N+1 -t N ) ^ ^ 

Let us now follow the time evolution of the atomic 
density operator in a time interval [0,t] for a sequence 
of atom detections: we assume that the first atom is de- 
tected at the position corresponding to <f>i and at time 
t\, the second atom at <f>2, t'l etc. and the 7V-th atom 
at <f>N, tjv- In the present case, starting from the initial 
state (4 ) and according to (8) and (9), this sequence of 
N detections prepares at time t the density operator 

?±l[\a + be i «*-»\*R atbt (<f>) (10) 

71 k = l 

Jo 27T 

We identify f { a P (<j>) oc Uk=i l« + &e^-^| 2 with the 
distribution of the relative phase of a state with decaying 
coherent amplitudes a t = ae~ Kt , b t = be~ Kt , and C is a 
time-independent normalization constant. 

The simplicity of Eq. (10) is due to the properties of 
coherent states: Firstly, they remain coherent but with 
a decaying amplitude under the time evolution with the 
effective Hamiltonian if e ff, i-e. e~ lH ^ fft \a)\b) oc |a t )|& t ). 
Secondly, they are eigenstates of the annihilation opera- 
tor and therefore also eigenstates of the jump superoper- 
ators: 

J iJ RabW = \a + be i ^-^\ 2 R ab (4>). (11) 

Note that f^\(i>) is independent of the jump times t^, 
which is a consequence of our assumption of equal damp- 
ing constants, i.e., equal transverse momenta of the two 
modes. 

The probability for observing the next count at posi- 
tion corresponding to the angle ip is thus given by 

PM" rp-W + beW-'WJH*). (12) 
Jo 27T 

Suppose now that that in a particular run an atom is 
counted at 4>N+i ■ The state immediately after the N+l— 
th detection is 

p t+dtX ' ^\a + be^^\ 2 fiPWR atbt (^, (13) 
Jo 27T 
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i.e., the relative phase distribution changes according to 

the map f^ +1) {<t>) |a + &e^ N + 1 -^| 2 /if ) W (except 
for a normalization constant). Initially, before the first 
atom is detected, the position (relative phase) distribu- 
tion is uniform and P(ip) = const. The (random) position 
of the first and subsequent detected atoms gradually de- 
termines the relative phase of the condensates, thereby 
breaking the initial symmetry of the distribution. 
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2. Evolution of the distribution of the relative phase 
as a function of the number of detected atoms 
N = 1, 2, . . . , 8 and N = 50 for two sample runs. Note, that 
the scaling of the y-axis is arbitrary. 



Figures 2 and 3 give two examples of trajectories which 
were obtained by simulating a sequence of counts assum- 
ing equal initial amplitude a = b (note that in that case 
the result is independent of the specific value of a = h) . 
In Fig. 2(a,b) we have plotted the evolution of the relative 

phase distribution fi^\<i>) after TV = 1, . . .,8 and after 
the 50-th detected atom. Typically, the phase distribu- 
tion is well-peaked after the first few counts. Asymptoti- 
cally the width of the distribution scales as A(f> oc 1/a/TV, 
and p t approaches a state with a well defined relative 
phase. Of course, the value of this phase varies randomly 



from run to run. If we average over many realizations, the 
interference pattern disappears, as expected from master 
equation (5). Figure 3 shows a simulation of the inter- 
ference pattern in the x—y plane of a single realization of 
the experiment for three different number of detections 
TV = 20,200 and 1000 (a,b,c). 
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FIG. 3. Buildup of the in- 

terference pattern for TV = 20, 200, 1000 atoms corresponding 
to (a), (b) and (c), respectively. 



Note that the state of relative phase R(<f>) is indeed 
a fixed point of the map defined by the simulation pro- 
cedure, since it is reproduced by any jump operator J^, 

Ji> Rab{4>) /Tr{J"^ R a b{4>)} = Rab{<t>), i- e -> R ab{4>) does not 
evolve apart from a decay of its amplitude. The distribu- 
tion of atoms in the screen corresponding R(<f>) is given 

by 

P(ip) oc \a + be^-^ \ = a 2 + b 2 + 2ah cos(V' - <t>) (14) 

and agrees perfectly with the interference pattern that 
one would expect intuitively from a state of fixed relative 
phase <f>. 

An interference experiment of the type described here 
prepares a state of the form (3) with a fixed relative 
phase between the two condensates. Note that this s- 
tate is diagonal in the total number of particles, i.e. 
does not violate superselection rules, and therefore it 
gives a consistent description of the experiment. Fur- 
thermore, the interference fringes for any possible out- 
come (<f>i, <f>2, ■ ■ ■ ,4>n) are completely consistent with the 
outcome produced if one had initially two coherent states 
\a) and \be^), with a phase difference <f> coinciding with 
the asymptotic (TV — > oo) phase difference of R(<f>). This 
means that for the description of any interference experi- 
ment with Bose condensates, one can simply assume that 
the initial state of each of the condensates is a coherent 
state with a random phase. Despite the fact that this 
state does not comply with the particle superselection 
rule, the results predicted in any experiment will be the 
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same as those given by a consistent theory such as the 
one presented here. 

The results obtained above can be also explained and 
generalized using the coherent state representation. In 
order to do that, we use the following property of Quan- 
tum Mechanics: if the state of a system can be written in 
a diagonal form with respect to a set of states |a 8 -) (that 
do not need to be orthogonal), i.e. p = Pj|aj)(aj|, 
then a single realization of a given experiment must yield 
the same result as if the state of the system is one of these 
states \ati). Obviously, the probability of obtaining the 
result corresponding to |a 8 ) is Pi. In our case, we can 
expand (a large class of) states that are diagonal in the 
Fock basis representation in terms of coherent states [14] 

p = J daP(a) j** ' ^\ ae ^)(ae^\, (15) 

and in particular the states (2). Consequently, a single 
realization of any experiment must give a result compat- 
ible with at least one of the particular coherent states 
\ae l ^ a ), which leads to the same conclusion as in the pre- 
ceding paragraph. This interpretation can also be used 
to show that two condenstates in Fock states may lead to 
interferences [8,9]. For any physical state of the conden- 
sate diagonal in the atomic number, we can write (15) or 
(1) (both are equivalent [15]). Thus, a single realization 
of an experiment has to be compatible with both Fock 
and coherent states. Since for two inital coherent states 
the probability distribution is of the form (14) (the co- 
herent states are fixed points of any jump operator 3^), 
for two Fock states it must have the same form. 

As a final comment, we note that the above discussion 
implies the existence of a phase standard for the relative 
phase of between Bose condensates. Let us assume that 
initially we have three condensates A, B and C which 
have been prepared in the product state po = p^ (dp^" 1 <X> 
p( c ) corresponding to an ensemble of coherent states | 
ae^") | be i( t">) | ce ! '^) with uniformly distributed random 
phases <f> a ,<f>b, and <f> c . First, observation of intereferences 
between A and B in a single run will fix the relative 
phase between A and B, a = <f> a — <f>b, and thus prepare 
the state | ae^") | be^"'") | ce 8 '^) with <f> a ,<f> c random. 
Second, observation of interference between B and C will 
prepare a state | ae^") \ be^"'") | ce 8 '^"' 3 ) with <f) a 
random. From this we can now predict the phase which 
will be observed in an interference experiment between 
A and C, 7 = a + [3, at least within the uncertainty of 
our knowledge of the phases a and [3. 

In summary, using a description based on continuous 
measurement theory we have shown how the dynamics 
in a single run of an interference experiment between t- 
wo Bose condensates prepares a state with relative fixed 
phase. Finally, we emphasize that our discussion explic- 
itly assumes that time scale of phase fluctuations due to 
dynamical properties of the Bose condensate [17] is much 
slower than the observation times necessary to build up 
an interference pattern. 



After this work was completed we have become aware 
of two related discussion of the phase of Bose condensate 
and its measurement by J. Dalibard and K. M0lmer [16]. 

This work was supported by the Austrian Science 
Foundation, the Marsden Fund contract GDN-501 un- 
der the auspices of the Royal Society of New Zealand, 
and the Deutsche Forschungsgemeinschaft. 



[1] M.H. Anderson et al. Science 269, 189 (1995). 

[2] C.C. Bradley, C.A. Sacket, JJ. Toilet and R. Hulet, Phys. 

Rev. Lett. 75, 1687 (1995). 
[3] K.B. Davis et al. Phys. Rev. Lett. 75, 3969 (1995). 
[4] For a review see Bose-Einstein Condensation, edited by 

A. Griffin, D.W. Snoke and S. Stringari (Cambridge U- 

niversity Press, 1995). 
[5] A.J. Leggett in Ref. [4] 

[6] An analogous situation exists in laser theory where the 
number of atomic excitations plus the number of photons 
is conserved in the rotating wave approximation. Thus 
for an intially empty cavity and a mixed state of atomic 
excitations the reduced density operator of the photon 
field will be diagonal in the photon number basis. 

[7] For are review and a complete list of references see P. 
Zoller and C. W. Gardiner in Quantum fluctuations, Les 
Houches, eds. E. Giacobino et al. (Elsevier, NY, in press). 

[8] J. Javanainen and S. M. Yoo, Phys. Rev. Lett. 76 , 161 
(1996). 

[9] M. Naraschewski, J. I. Cirac, H. Wallis, A. Schenzle, and 
P. Zoller, Phys. Rev. A, (to be published). 

[10] For simplicity we ignore the effects of gravity. 

[11] This assumption is not essential. Employing a P- 
representation Eq. (2), the following derivations are read- 
ily generalized to arbitrary initial states. 

[12] This becomes apparent if we rewrite (3) as 

N=0 

where 

N / \ ~ 

= (a 2 +6 2 )-f J2 ( I, ) « {N - n) t> n e-^\(N-n)) a \n) b , 

n = ^ ' 

is a state with N particles. 
[13] M. Naraschewski, diploma thesis, Universitat Miinchen 
(1993). 

[14] C. W. Gardiner, Quantum Noise (Springer- Verlag, 
Berlin, 1991). 

[15] The coherent state is an overcomplete basis, which allows 
to write a given density operator in a diagoanal form in 
two different sets, the coherent states and the Fock states. 

[16] Talks given by J. Dalibard and K. Mollmer, Confer- 
ence on Collective Effects in Ultrcold Atomic Gases, Les 
Houches (April 1 - 5, 1996). 

[17] M. Lewenstein, private communications. 



4 



